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Abstract 

We discuss nonanomalous R-symmetry in the supersymmetric grand unified the- 
ories. In particular, we explore anomaly-free solutions predicting the gravitino 
mass in the range of 10~^ eV^ 777,3/2 ~1 TeV when the ;U-parameter is fixed to 
be ^ ~ 1 TeV. In the minimal SU(5) GUT, we have shown that ;U ~ 1 TeV is 
obtained only if the gravitino is ultralight with mass 771.3/2 ~ 10^3 eV. If extra 
fields 5 © 5* or 10 © 10* are introduced, many solutions predicting 7773/2 ^ 10"'^ eV 
are found. The R-parity is violated due to the vacuum expectation value of the 
superpotential, but it is controlled by the discrete R-symmetry. We find that the 
R-parity violating couplings are naturally suppressed much below the experimental 
bounds for some charge assignments. These charge assignments predict light grav- 
itino with masses of order 0(10"^ eV)-C'(l MeV). These discrete R-symmetries can 
be considered as solutions to the /7-problem in low energy supersymmetry breaking 
models such as the gauge mediation. 



1 Introduction 



A discrete R-symmetry Z^r often appears as a remnant of the rotational symmetry of 
the compactified extra space in higher dimensional supergravity or string theory [1, 2]. 
This discrete R-symmetry should be nonanomalous since this is a gauge symmetry. An 
R-symmetry plays a crucial role in the phenomenology of supersymmetric (SUSY) theory. 
First, it can suppress the cosmological constant compared to the Planck scale. Second, 
the SUSY-invariant mass term (the /i-term) of the Higgs chiral multiplet can be forbidden 
so that the Higgs mass not be the Planck scale. If an R-symmetry breaking is related 
to SUSY breaking, the Higgs chiral multiplet can obtain a mass of the order of the 
gravitino mass m^i2 — 1 TeV by the Giudicc-Masiero (GM) mechanism [3]. Third, an R- 
parity forbids the dimension-four baryon and lepton number violating operators causing 
too rapid proton decay [4, 5]. These observations motivate us to ask whether wc can 
find a nonanomalous discrete R-symmetry with the above properties. In a paper by 
Kurosawa, Maru and Yanagida [6], nonanomalous discrete R-symmetries in the minimal 
SUSY standard model (MSSM) and the SUSY grand unified theory (GUT) were found 
under the situation that the GM mechanism works. These solutions can also forbid the 
dimension-five baryon- and lepton-number violating operators. Furthermore, extra fields 
5 ® 5* with the mass of order 1 TeV, which can be testable in collider experiments, are 
predicted from the anomaly cancellations in the GUT case. 

In this paper, we consider nonanomalous R-symmetries without the GM mechanism.^ 
In particular, we consider that the //-term is induced by the vacuum expectation value 
(VEV) of the superpotential (W), assuming that a fractional power of (1^) is allowed in 
the superpotential. The fractional power makes it possible to obtain the correct size of 
the /x-term even for gravitino mass smaller than the electroweak scale. We find that the 
/i ~ 1 TeV is obtained only if the gravitino is extremely light as m3/2 — 10~^ eV in the 
minimal SU(5) GUT. If extra fields such as 5 © 5* or 10 © 10* are added to the minimal 
SU(5) GUT, we find many charge assignments predicting 10~^ eV. These charge 

assignments can be considered as solutions to the /i-problem in low energy SUSY breaking 
models such as the gauge mediated SUSY breaking (GMSB) models [8].^ 

In our framework, the R-parity is in general violated due to the fractional powers of 
(W). However, the R-parity violating couplings are well controlled by the symmetry. In 
fact, it turns out that R-parity violation is small enough for some charge assignments. It 
is further shown that the dimension-five baryon- and lepton-number violating operators 
are naturally suppressed. In these charge assignments the gravitino masses are predicted 

^We do not consider a possibility of anomaly cancellations by Green-Schwarz mechanism [7]. 

^Solutions to the ^-problem in the GMSB models have been reported so far in Rcfs. [9, 10, 11, 12]. 
Also interesting mechanisms for the /x-problem in various mediation mechanisms of SUSY breaking are 
recently proposed [13, 14]. 
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in the range m3/2 ~ (9(10^'^ eV)-(9(l MeV), and hence the gravitino is the hghtest 
supersymmetric particle (LSP). Its Ufetime is much longer than the age of the universe, 
and the gravitino can be the dominant component of the dark matter. 

This paper is organized as follows. In the next section, after discussing the general 
constraint on the power of the superpotential, we search for anomaly-free solutions of 
the minimal SU(5), and the minimal SU(5) with 5 ® 5* or 10 © 10*. In our analysis, 
constraints on the R-parity violating operators from the proton decay and neutrino masses 
are taken into account. A brief discussion on the cosmology of the gravitino LSP with 
R-parity violation is also given. The last section contains a summary of our paper. In 
the appendix, a subtle issue between the fractional power and the discrete symmetry is 
discussed. 



2 Discrete R-symmetry in GUT 

In this section, we consider the discrete R-symmetry in the GUT. Before discussing 
anomaly cancellations in detail, we give a general constraint on powers of the VEV of 
the superpotential in the next subsection. 

2.1 General constraint on the powers of W 

We assume that the //-term is generated from 

where y is a non- negative number and Mp = 2.4 x 10^^ GeV is the reduced Planck scale. 
Here and hereafter, we omit coefficients of order unity. Then, the following conditions 
should be satisfied, 

2ay + h + h = 2a (mod A^) , (2) 
h + h ^ 2a (mod N) , (3) 

where h, h and a denote the R-charge of H, H and the Grassmann coordinate 9, respec- 
tively and they are all integers. (See Table 1.) The second condition (3) is necessary to 
forbid the Higgs mass term with Planck scale. Prom Eq. (1), we can predict the gravitino 
mass because 



and hence 

^3/2 



Mr 



Mp ~ GeV , (5) 
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T F N H H e 


SU{5) 
Znr 


10 5* 1 5 5* 

t f n h h a 



Table 1: The matter content of GUT. Z^r charges of the fields denote those of the scalar 
components and are integers. We take the R-charge of the Grassmann coordinate 9 to be 
an integer a. 



where /j, ~ TeV is assumed throughout this paper. Eq. (5) can be rewritten as 

log(Mp/^) 



(6) 



log(Mp/m3/2) ' 

Prom this equation, we can derive lower and upper bounds on y. We assume 

1713/2 < pi ^Mp. (7) 

Thus, Eqs.(6) and (7) lead to 

0<y<l. (8) 

Here, one might wonder if the fractional power y < 1 is incompatible with the discrete 
symmetry. This issue is briefly discussed in the Appendix. As can be seen from Eq. (4), 
the fractional power y < 1 is crucial to obtain the correct size of /x-term for a gravitino 
mass smaller than the weak scale. 

Besides the constraint in Eq.(8), there is a lower bound on y coming from the lower 
bound on the gravitino mass m3/2- In the GMSB model, when we fix the soft mass 
scale msoft, the SUSY breaking F-term is bounded from below as -\/F^O(10 TeV) x 
(msoft/100 GeV) in order to avoid the negative mass squared for the scalar field in the mes- 
senger sector [8]. This leads to a lower bound on the gravitino mass, 777.3/2 ^ 0(0.01 eV). 
Even fighter gravitino mass 771.3/2 ~ O{10~^ eV) can be allowed in some SUSY breaking 
models in higher dimensional spacetime [15]. For 777.3/2 ~ 10~^ eV, one can see from 
Eq. (6) that y can be as small as y ~ 1/2. Prom the above arguments, we impose the 
following constraints on the parameter y: 

\<y<^- (9) 

2.2 Anomaly cancellation 

We are now at the position to discuss the anomaly cancellation. Let us first take the 
minimal SU(5) GUT. Its matter content is described in Table 1. Z^r charge of the 
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fields, which is taken to be generation independent for simplicity, denotes those of the 
scalar component and are integers. Note that we take the R-charge of the Grassmann 
coordinate 9 to be an arbitrary integer a. The Yukawa coupUngs and the Majorana mass 
term are given by^ 

W = TTH + TFii + FNH + -MmN'^ , (10) 

where is a Majorana mass for the right-handed neutrinos. For Eq. (10) to be allowed, 
the corresponding R-charges have to satisfy 

2t + h = 2a (mod N) , (11) 

t + f + h = 2a (mod A^) , (12) 

f + n + h = 2a (mod N) , (13) 

2n = 2a (mod N) . (14) 

Anomaly cancellation conditions for ZNji[SU{3)c]'^ and ZNji[SU{2)i]'^ are [16] 

Q _ AT 

ZNn[SU(3)c]'' - -{3it-a) + if-a)} + 3a^-k (k e Z) , (15) 
Znr[SU{2)l]^ = ^{9{t-a) + 3{f-a)} + ^{{h-a) + {h-a)} + 2a 

= ^k' (k'eZ). (16) 
These conditions are simplified to 

h + h = 2a (mod — ) , (17) 

N 

h + h ^ a (mod — ) . (18) 

Eqs. (17) and (18) lead to 

h + h = ^a (mod iV) . (19) 

Substituting this back into (17) or (18) results in 6q; = (mod A^), or equivalently 

6a = Nk (A; e Z) . (20) 

Taking into account 0<q;<A^, 0<A;<6is obtained. If we take k — 3^ then 2a — N, 
and 

h + h^4:a^2a (mod 2a) , (21) 



We have suppressed 0{1) coefficients of the terms. 
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is derived and contradicts the condition (3). Thus, we find k = 1, 2, 4 and 5. In other 
words, N is classified as, 

N = 6a, 3a, 3 (|) , 6 (|) . (22) 

First, N = 6a case is considered. From Eqs. (2) and (19), 

2ay + Aa^2a (mod 6a) ^ y ^ -1 + 3n (n e Z) . (23) 

This has no solution for 1/2 < y < 1. Second case is N — 3a. Prom Eqs. (2) and (19), 
we obtain 

3 

2ay + 4:a^2a (mod 3q;) ^ y = -1 + -n (n e Z) . (24) 

This has a solution y = 1/2 for n = 1. Third case is = 3 (q;/2) = 3q;' {a — 2a', a' e Z). 
From Eqs. (2) and (19), we obtain 

2ay + Aa = 2a (mod 3a') , (25) 

3 

<^ Aa'y + 8a' = Aa' (mod 3a') ^ y ^ -1 + -n (n e Z) . (26) 

This also has a solution y = 1/2 for n = 2. The last case is = 6 (a/5) = 6a' (a = 
5a', a' e Z). From Eqs. (2) and (19), we obtain 

2ay + 4a = 2a (mod 6a') , (27) 

3 

<^ lOa'y + 20a' = 10a' (mod 6a') ^ y ^ -1 + -n (n e Z) . (28) 

5 

Taking into account 1/2 < y < 1, 5/2 < n < 10/3 is obtained. This has a solution n = 3, 
but it does not correspond to the minimum solution for positive y. The dominant /x-term 
comes from the term with y = 1/5 (i.e. n — 2), and it would cause // » 1 TeV unless 
"7.3/2 <^ 10~^ eV. Thus, this case has no solution. 

Therefore, the correct size of the //-term can be obtained for y = 1/2 when N — 3a and 
N — 3{a/2), which predict an extremely hght gravitino with mass 777,3/2 — O{10~^ eV). 

Bciryon- and lepton-number violating operators 

We have shown that the correct size of //-term can be obtained for ultralight gravitino 
7773/2 ^ 10'^ eV, by means of the fractional power of the superpotential's VEV, (W)^^^. 
However, if we allow general interaction terms including the fractional power of (W), 
there appear baryon- and lepton-number violating operators as well. Therefore we have 
next to consider constraints on these operators. 
Let us first consider the following superpotential 

^"^{m) ^^^^ 
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which includes the so-called bilinear R-parity violation, W 
given by 



\Mp) 



p ■ 



fiiLiHu. The coupling is 



(30) 



The bilinear R-parity violation can generate the neutrino mass [17] which explains the 
atmospheric neutrino oscillation if /i/yU ~ (9(10~^-10~^) [18, 19]. In other words, p, should 
be smaller than C(10^^ GeV) in order to avoid a too large neutrino mass. 

The power z is determined by the symmetry: 



2az + f + h = 2a (mod N) , 

2az = n (mod N) , 
n + rN . ^, 
2a 



(31) 



where we have used Eq. (13). From Eq. (14), the charge of the right-handed neutrino 
should be either n = a or n = a + N/2 (mod N). li n — a, the minimum non- negative 
z is given by z = 1/2 since N > a. In this case, the bilinear coupling is given by 



Mp 



1/2 



Mp ~ 10^ GeV 



^3/2 

10-3 eV 



1/2 



(32) 



This generates too large neutrino mass and hence is excluded. Thus, the charge of the 
right-handed neutrino should be ri = a -|- N/2 (mod A^). This also means that the case 
oi N — odd is excluded. UN— even and n — a + N/2 (mod N), Eq. (31) gives rise to 



5 + 6r 
4 

7 + 6r 



for N = 3a , 



for 7V = ^. 

2 



(33) 



The minimum non-negative z are given by z = 5/4 (r = 0) and z = 1/8 (r = — 1), 
respectively. Therefore, the bilinear R-parity violating coupling is given by 



/^3/2 



Y^Mp 3x 10-20 GeV (j^pf^)^^^ N^3a, 



^3/2 

Mp 



1/8 



Mp ~ 4 X 10^^ GeV ( ^T^X ^^ for N 

\w~'^ eV/ 



(34) 



3a; 
T 



Thus, the = 3(q;/2) case is clearly excluded. On the other hand, for N — 3a the con- 
tribution to the neutrino mass from the R-parity violation is extremely small. Therefore, 
the dominant contribution to the neutrino masses are understood to be generated by the 
standard seesaw mechanism [20] . 
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Next we consider the trilinear R-parity violation caused by the following superpotential 

which includes both of the baryon- and lepton-number violating operators, XijkLiLjE^, 
X'^j^LiQjDl and X'lj,^U^DjDl. The power z' is determined by 

2az' + t + 2f^2a (mod N) , 

2az' ^2a + n (mod N) , (36) 

where we have used Eqs. (12), (13) and h + h = 4a. For N = 3a {= even) and n = 
a + N/2 (mod N), z' is given by 

q _|_ 5^' 

= ^ (r'eZ). (37) 

The minimum non- negative z' is given by 2;' = 3/4 {r' — —1) and the effective coupling 
becomes 



\ 3/4 

) , (38) 



which easily satisfy the constraint from the proton decay, X'ujX'uj $2x10 ^^(mgoft/lOO GeV)^ 
(J = 2,3)[21]. 

Let us also discuss the dimension five operator W ~ {l/Mes)TTTF, which causes 
the proton decay [5] . In order to suppress the proton decay rate below the experimental 
bound, the effective mass scale should be M^s > 0(10^^ GeV) [22], i.e., much larger than 
Mp. Notice that the R-parity cannot forbid this operator. In our framework, the operator 
is given by 

W c^ — TTTF ~ — f TTTF = -^TTTF , (39) 

and z" should satisfy 

2az" + 3t + f^2a (mod N) , 
2az" = 2a (mod N) , 

2a 2 ^ ^ ^ 

where we have used Eqs. (11), (12), h + h — Aa, and N — 3a. The minimum non-negative 
z" is given by z" — 1 (r" = 0), and hence the effective mass scale Mes ~ Mp{Mp/ms/2) ~ 
10^^ GeV(10~^ eV/ms/2) is much above the experimental bound. Therefore, the discrete 
R-symmetry naturally suppresses the dimension five proton decay operator. 
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t 


/ 


n 


h 


h 9 


p 


1 


1 


1 





a 


V 


1 


-3 


5 


-2 


2 


A 





-1 


1 





1 



Table 2: The charges of the generators P,V,A. 



explicit Znr charge assignments 

Before closing this subsection, we comment on the explicit form of Znr. Nonanomalous 
R-symmetries are represented by Z^r — P°'V^A'^ [a, (3,^ e Z) [16] where the generators 
P, V and A are summarized in Table 2. One can easily check that the charge assignments 
in Table 2 satisfy Yukawa conditions (11)-(13).'^ For N — ?>a {— even) and n — a -\- 
N/2 (mod A^), we obtain 

N iV 3 

n = Q; + 5/3 + 7 = Q;H (mod N) ^ 5/? + 7 = — = -a (mod N) , (41) 

while Eq. (19) leads to 

-2(5 + (2/? + 7) = 7 = 4q; (mod N) . (42) 

From Eqs. (41) and (42), we obtain 

Z-io.R = {PA^'yvP, (43) 
5 

5/3 = "2*^ (mod N — ?,a) 
5 

= --a + 'iam {m E Z) . (44) 

Therefore, 

a : /3 : = 10 : (-5 + 6m) : 30 , (45) 

and the explicit forms of the Z^r are given by 

Z3o« = (PA^)iV6— 5 (m=l,2,3,4), (46) 
Z^R = (PA'l)V^ (47) 



^By imposing 3 conditions (11)-(13) on the 6 parameters t, f, n, h, h and 9, the charge assignments 
can be represented in terms of 3 parameters a, (3 and 7. 
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2.3 Introducing 5 © 5* 

In this subsection and the next, we consider the possibihty that 10~^ eV ^ m3/2 ^ 1 TeV 
is predicted by introducing a pair of 5 © 5* or 10 © 10*. If we introduce a pair of 5 © 5* 
with Zjvj? charges ^ and ^, respectively, the anomaly cancellation conditions are modified 
as 

Znr[SU{3)c]' = ^{3{t-a) + {f-a)} + ^{{^-a) + {^-a)} + 3a^^k {k E Z) 
44> 3{h + h)-{^ + 0=^a (mod N) , (48) 
Znr[SU{2)l]^ = ^{9{t-a) + 3{f-a)} + ^{h + h-2a) 

1 - N 

^ 2{h + h) = ^ + ^ (mod A^) . (49) 
Prom these conditions, 

h + h^Aa (mod N) , (50) 

^ + ^^Sa (mod N) , (51) 

are obtained. At this stage, is undetermined. In order to fix A^, let us take into account 
the mixed gravitational anomaly cancellation [16], 

ZjvR [gravity]^ = 30{t - a) + 15(/ - a) + 3{n - a) + 2{h - a) + 2{h - a) 
+5(^ - a) + 5(^ - a) + (8 + 3 + l)a - 21a 
= -23a (mod A^ or N/2) . (52) 

Mod N or N/2 depends on whether N is odd or even. In the following, we consider 
whether the mixed anomaly can be canceled without introducing singlets. 
If A^ is odd, 

23 

N^—a {k^ 1,2, ■■■,22). (53) 



k 



Eq. (5) tells us 



because 



18.4- 



15.4 



m3/2 ~ 10"" -1+1" GeV (n e Z) , (54) 



23 

2ay + 4:a = 2a {mod 23a/k) y = -1 + —n . (55) 

Zrv 

Taking into account that 1/2 < y < 1 and y is a minimum of non- negative value, we 
obtain 

A; = 6, 7 (for n = 1) , (56) 
k = 12, 13, 14, 15 (for n = 2) . (57) 
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N : odd 












(n 


'"'■if 2 


( n fc) 


'"'■if Z 


(n fc) 


'"'6/2 


(1,6) 


40 GeV 


(2, 12) 


40 GeV 


(2, 14) 


2.9 keV 


(1,7) 


2.9 keV 


(2, 13) 


24 MeV 


(2, 15) 


0.035 eV 


N : even 












(1, 12) 


40 GeV 


(2, 24) 


40 GeV 


(2, 28) 


2.9 keV 


(1, 13) 


24 MeV 


(2, 25) 


1.2 GeV 


(2, 29) 


14 eV 


(1, 14) 


2.9 keV 


(2, 26) 


24 MeV 


(2, 30) 


0.035 eV 


(1, 15) 


0.035 eV 


(2, 27) 


330 keV 







Table 3: The gravitino mass for GUT with 5 ® 5*. 



The gravitino mass is summarized in Table 3. 
If N is even, 



46 

N^—a (fc = l,2,---,45), (58) 
k 

23 

2Q;y + 4Q; = 2Q; (mod 46q;/A;) ^ y = -1 + — n , (59) 

k 



and Eq. (5) leads to 



18.4- 



15.4 
"23- 



ms/2 ~ 10 -i+t" GeV (n e Z) . (60) 



Taking into account that 1/2 < y < 1 and y is a minimum of non- negative value, we 
obtain 

A; = 12 ~ 15 (for n = 1) , (61) 
/c = 24 ~ 30 (for n = 2) . (62) 

The gravitino mass in this case is also summarized in Table 3. 

Now let us turn to discuss the baryon- and lepton-number violating operators. First, 
the bilinear R-parity violating coupling ji is given by Eq. (30), with a fractional power z 
in Eq. (31). One can show in the same way as before that the case of n = a is excluded 
since it results in z = 1/2, which induces too large neutrino masses. UN — even and 
n^a + N/2 (mod N), Eq. (31) and Eq. (58) give rise to 

fe + 23 + 46r 

z = — (r e Z) . (63) 

The minimum non-negative z is given by r = for k — 12-15, and r — —1 for k — 24-30. 
The resultant bihnear couphngs are given by /2 ~ (8 x 10~^-7 x 10~^^) GeV for k — 12-15, 
and fl ~ (10^^-10^^) GeV for k = 24-30. Hence, k = 24-30 cases are excluded. Namely, 
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k 


m3/2 






14 
15 


2.9 kcV 
0.035 eV 


6 X 10-14 GcV 

7 X 10-19 GeV 


6 X 10-1^ 

7 X 10-22 



Tabic 4: The gravitino mass for GUT with 5 © 5* and R-parity violating couplings p, and 
Acfr- N = even and n = a + N/2 (mod N). Only the cases which predict sufficiently 
small R-parity violation are listed. 

the remaining charge assignments are N — {46/ k)a — even, n — a + N/2 (mod N), and 
k = 12-15. 

Next we consider the trilinear R-parity violating coupling W ~ XesTFF. The effective 
couphng Aeflf is given by Aeff ~ (m3/2/Mp)^'. The power z' is again determined by Eq. (36), 
2az' — 2a + n (mod N). (Notice that h + h — Aa is satisfied also in the present case. See 
Eq. (50).) Therefore, for N = {A6/k)a = even and a + N/2 (mod N), z' is given by 

, 3/c + 23 + 46r' , , - 
z'^ ^ (r'eZ). (64) 

The minimum non-negative z' is given by r' = for k = 12-15, and the trilinear couplings 
are Acg ^ 8 x lO^i^ {k = 12), 5 x lO^^^ {k = 13), 6 x lO^i^ (A; = 14), and 7 x 10^22 
{k = 15). Thus, the cases of A; = 12, 13 are excluded by the constraint from the proton 
decay, X[,j\'l,^ < 2 x 10-2^ (msoft/100 GeV)2 {j = 2, 3) [21]. 

The R-parity violating couplings are listed in Table. 4 for the remaining charge as- 
signments k — 14, 15 together with the gravitino masses. We find that the neutrino 
mass induced by these R-parity violations is too small to explain the mass scale observed 
in the atmospheric and solar neutrino oscillation experiments, and hence the dominant 
contribution to the neutrino masses should be generated by the seesaw mechanism. 

As for the dimension five operator W ~ {l/Mes)TTTF, the effective mass scale is 
given by Eq. (39), M^s ~ Mp{Mp/ms/2y" ■ It is easy to show that the power is again 
given by z" — 1 for = {46/ k)a = even, n = a + N/2 (mod A^) and k = 14, 15. Thus, 
the effective mass scale M^s ~ Mp{A{p/ 111^12) ~ 10'^^ GeV(l keV/m3/2) is much above 
the experimental bound and the proton decay via the dimension five operator is naturally 
suppressed by the discrete R-symmetry also in these charge assignments. 

Finally, we comment on the explicit form of Z^p = P'^V^A"' {a, f3,'~f E Z). (See Table 
2). We only consider the cases of A^ = even, n = a + N/2 (mod A^) and k = 14, 15, 
since other cases are excluded by the constraints on R-parity violation, as we have shown. 
Thus, 

A^ A^ 

fl^a + 5(3 + -f^a-{ (mod A^) <^ 5/3-^7 = — (mod A^) . (65) 

2 2 
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Eq. (50) leads to 

-2(3 + (2/3 + 7) = 7 = 4q; (mod N) . (66) 
Prom Eqs. (65) and (66), 



N 

5/3 + 4a = — (mod N) , (67) 



is obtained. Thus, we have 



^faii = {PA'rv''. (68) 
.n. 23 , , 46 , 23 46 

5/3 + 4q; = —a (mod — a) = —a + —am (m G Z) . (69) 
k k k k 

Therefore, 

a : /3 : iV = 5A; : (46m + 23 - Ak) : 230 , (70) 



and we obtain 



Z230R = (PA^)™\/46™-33 ^j.. = 14 , ^ = 1, 2, 4, 5) , 

(PA^)^5l^^^™-3^ (A; = 15, m = 1,3,4,5). (71) 
ZmR = (P^^)'V2^ {k^U), {PAyv^^ (A; = 15). (72) 



2.4 Introducing 10 10* 

Next, we consider the case with ^(10) ® ^(10*). The anomaly cancellation conditions are 
modified as 

3 - 3-3 TV 

Znr[SU{3)c? : -^{h + h) + 3a+-{^ + + ^Nk" - 3a ^ -k , (73) 

3 N 

Zmr[SU{2)l]' : -{h + h) + a + -{^ + 0-^Oi^^k'. (74) 
These are simphfied to 

3{h + h) = 3(^ + (modiV), (75) 
2{h + h) = 3(C + - 4a (mod TV) , (76) 

and then 

h + h^4a (mod N) , (77) 
3(^ + ^) = 12q; (mod N) , (78) 
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are obtained. In order to fix A^, let us take into account the mixed gravitational anomaly 
cancellation, 



Z^ij [gravity]^ = 30{t - a) + 15(/ - a) + 3{n - a) + 2{h - a) + 2{h - a) + 10(^ - a) 
+10(^ - a) + (8 + 3 + 1)q; - 2la 
= + - 37a (mod N or N/2) . (79) 

In the same manner as the ^(5) ® case, one can show that the n ~ a case is 

excluded because it would generate too large neutrino mass from the bilinear R-parity 
violation. Therefore, we consider the case N — even and n — a + N/2. Then, Eqs. (78) 
and (79) lead to 

3ZjvR [gravity]^ = -99q; (mod N/2) , (80) 
198 

^N = —a (A;= 1,2,---,197). (81) 



Prom Eq. (2), we obtain 



99 

y^-l + —n (neZ). (82) 

fV 



Taking into account that 1/2 < y < 1 and y is a minimum of positive values, 

= 50 ~ 66 (for n=l), (83) 
A; = 100 ~ 132 (for n = 2) , (84) 

are obtained. 

Some of these charge assignments cause too large baryon- or lepton-number violation 
via the R-parity violation. The orders of magnitudes of bilinear (/i) and trilinear (Aes) 
couplings can be estimated in the same way as 5 ® 5* case. In Table. 5, we show the 
gravitino mass and these couplings for the cases in which the R-parity violations are 
below the experimental bounds {k = 57-66). In these cases, the dimension five operator 
W ~ {l/M^s)TTTF are naturally suppressed as M^s > 10^^ GeV like the 5 © 5* case. 

Finally, we comment on the concrete form of Z]sir. From Eqs. (77), (81) and n = 
a + N/2, we obtain 



{PAyvf^ , (85) 

99 198 99 198 

5/9 + 4q; = —a (mod — ~;~'^) = ~r(^ "I — —Oim (m e Z) , (86) 
k k k k 

and hence 

a: P:N^5k: (198m + 99 - 4A;) : 990 . (87) 
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k 


'''■3/2 




Aeff 


O ( 


O.Z iVlt: V 


i X 


1 n 
iU 


-11 




( X 


1 n- 
iU 


-14 


Oo 


/ion LroA^ 


y X 


iU 


-12 


LrC V 


y X 


iU 


-15 


oy 


AO 


1 X 


1 


-12 


vjrc V 


1 X 


1 0" 


-15 


60 


5.2 keV 


1 X 


10 


-13 


GeV 


1 X 


10- 


-16 


61 


490 eV 


1 X 


10 


-14 


GeV 


1 X 


10" 


-17 


62 


41 eV 


8 X 


10 


-16 


GeV 


8 X 


10- 


-19 


63 


2.9 eV 


6 X 


10 


-17 


GeV 


6 X 


10- 


-20 


64 


0.18 eV 


4 X 


10 


-18 


GeV 


4 X 


10- 


-21 


65 


0.95x10-2 eV 


2 X 


10 


-19 


GeV 


2 X 


10- 


-22 


66 


0.4x10-=^ eV 


9 X 


10 


-21 


GeV 


9 X 


10" 


-24 



Table 5: The gravitino mass for GUT with 10 ® 10* and R-parity violating couplings p, 
and Aefif. N — even and n — a + N/2 (mod N). Only the cases which predict sufficiently 
small R-parity violation are listed. 

For k = 66, it reduces to Z^qh or Zqji given in Eqs. (46) and (47). For k = 57-65, we 
obtain 

-- 63) , (88) 

= 57), (pA^)100y66m-47 (^^ 60), (89) 

(other k) , (90) 

where m — 1,2,3,4,5. Among these cases, for 198m + 99 — 4A; = 5m' (m' e Z), they 
reduce to Z22R, Zqqr, and ^iqsr- 

2.5 Gravitino LSP with R-parity violation 

As can be seen from Table. 4 and 5, the predicted masses of the gravitino are 0(10"^ eV- 
1 MeV). Therefore, the gravitino is the LSP. Here, let us comment on the cosmology 
of this gravitino LSP. Gravitino LSP dark matter without R-parity was investigated in 
Ref. [23] under the assumption that the R-parity violation is the dominant contribution to 
the neutrino masses. According to them, we have found that the lifetime of the gravitino 
is much longer than the age of the universe,^ since the decay rate is suppressed by the 
small R-parity violating coupling in addition to the Planck scale. (Notice that the R- 
parity violation in our scenario is even smaller than that considered in Ref. [23]. Thus, 
the lifetime of the gravitino in our case is much longer than that in their case.) For such 
a long lifetime, the flux of the diffuse gamma ray generated by the gravitino decay is 

^Thc lifetime of the LSP gravitino could be shorter than the age of the universe [24] if there is a large 
trilinear coupling A ~ 0(0.1-1) close to the experimental bound and if the gravitino mass is relatively 
large, 1713/2 ^ 0{1 GeV). However, this is not the case in our scenario. 



ZnoR = (PA'f'v'''--'' {k 

y /' P /1 4 ^ S/c T / 198m+99-4fc 
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smaller than the observed value [23]. Therefore, for m3/2 > 0{1 keV), the gravitino can 
be the dominant component of the dark matter in spite of the presence of the R-parity 
violation.^ 

As for the next-to-lightest SUSY particle (NLSP), it can decay either (i) into gravitino 
via the usual R-parity conserving interaction, or (ii) into the standard model particles via 
R-parity violating couplings. We have found that the partial decay rate of the latter one 
is much smaller than that of the former one, since the R-parity violations are extremely 
small. (See Table. 4 and 5.) Thus, the NLSP decay rate is determined by the former 
channel. If the decay occurs during or after the big-bang nucleosynthesis (BBN), i ~ 1- 
100 sec, it might spoil the success of the BBN [25]. However, for 1113/2 < 0{1 MeV) the 
lifetime of the NLSP is shorter than 1 sec and this problem is avoided. 

3 Summary 

In this paper, we have studied nonanomalous discrete R-symmetry in GUT without im- 
posing the Giudice-Masiero condition. In the minimal SU(5) GUT, ~ 1 TeV is obtained 
only if the gravitino is ultralight as 1713/2 — 10^^ eV, and we find simple solutions Zg^ 
and Z^Qpi^. If a pair of 5 © 5* or 10 © 10* are added to the minimal SU(5) GUT, we can 
find many solutions predicting 112^/2 ^ 10^"^ eV. Here, we comment on the mass of these 
additional multiplets, ^ + C,, which can be estimated since their charges are determined 
by the anomaly cancellation condition. [See Eqs.(51) and (78).] The effective operator 
which induces the mass of ^ + ^ is given by W ~ {{W) /M%)y'Mpii. We have checked 
that the mass ~ [1x13/2/ MpY Mp is larger than the electroweak scale in all cases we 
have discussed. 

Since the fractional power of {W) is considered in this paper, R-parity is necessarily 
violated, but R-parity violating couplings are controlled by the symmetry. In fact, the 
couplings were found to be small enough to avoid the constraints from proton decay and 
neutrino masses for some charge assignments. Furthermore, it has also been shown that 
dimension-five baryon- and lepton-number violating operators are naturally suppressed. 
Therefore, the proton stability is ensured by the symmetry. 

Low energy baryon- and/or lepton-number violating interactions might cause a dif- 
ficulty for baryogenesis since they might wash out the baryon asymmetry together with 
the sphaleron [26] process. However, we found that this is not the case for our scenario 
because the R-parity violating couplings are so small that their interactions have never 
been in thermal equilibrium. 

The predicted gravitino masses were found to be in the range 1113/2 ~ (9(10~^ eV)- 

^For m^/2 < 0{1 keV) the gravitino cannot be a cold dark matter, and hence the dominant component 
of the dark matter should be another particle or object. 
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0{1 MeV) and the gravitino is the LSP. Since the hfetime of the gravitino is much longer 
than the age of the universe, the gravitino can be the dominant component of the dark 
matter. 

As for the neutrino mass, there are two contributions in our framework, i.e., the seesaw 
mechanism and the R-parity violation. The R-parity violation can explain the neutrino 
mass scale if fi/fx ~ 0(10-^-10-^) and/or Aeff ~ 0(10-^-10-^) [19, 27, 28]. However, 
the predicted values of these couplings are either much larger or much smaller than these 
ranges. The charge assignments predicting too large couplings are excluded. Therefore, 
the neutrino mass scale should be generated by the seesaw mechanism, and contributions 
from the R-parity violation gives only tiny perturbation to it. 

The discrete R-symmetry considered in this paper can explain the order of magnitude 
of the //-term for hght gravitinos. Hence, they can be considered as solutions to the 
//-problem in low energy SUSY breaking models such as the gauge mediation. 

Note Added 

1 In our scenario, the predicted masses of the gravitino are very small. We would like to 
stress that, if the gravitino is hghter than about 100 keV, there is an interesting possibihty 
to detect the slow decay of the lightest neutrahno into gravitino in future coUider [29] . 

2 If the gravitino is lighter than the proton, proton can decay into the gravitino via the 
R-parity violating coupling \"U'^D'^D'^. This leads to a stringent limit on the couphng 
as A" < 10"^^(m3/2/eV) [30]. We found, however, that all the charge assignments which 
satisfy the constraints discussed in the text also satisfy this constraint, and hence the 
conclusion does not change. (We are grateful to Ryuichiro Kitano for pointing out this 
constraint.) 

3 5-parameter in our case is at most of order loop suppression factor times gaugino 
mass, which is induced at 2-loops in gauge mediation. This implies that large tan/3 is 
preferred. (We wish to thank Stephan Huber for useful comments and discussion.) 
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Appendix 



In this appendix we briefly comment on the fractional power and discrete symmetry. 
Although we have derived Eq.(2) from Eq.(l), it is nontrivial in the case of fractional 
power y < 1 since the R-charges are defined under mod N in the framework of discrete 
Znr symmetry. 

Let us first discuss the vacuum expectation value of the superpotential, (W) . Suppose 
that a gauge singlet operator X has a Z^r charge 2a and it develops a vacuum expectation 
value {X) <^ Mp. Under the Z^r symmetry, it is also possible to consider that X has a 
charge 2a + nN with n e Z. Then, in general, the superpotential can have the following 
vacuum expectation value: 



/ (X\\ Sa+nJV 

w = E Ml, (91) 

n>-2a/N J 

where we have renormahzed the mass dimension of X operator to be 1. If all of the coef- 
ficients c„ are of order one, the vacuum expectation value of the superpotential becomes 
Planck scale, which is inconsistent with almost vanishing cosmological constant or low en- 
ergy SUSY. Thus we expect that the right-hand side of the above equation is dominated 
by a certain term with n — Uq: 

2a 

/ (X) \ 2a+no^ 

(W) ^ Cno(^^j M|. (92) 

The number no is likely to depend on the operator X as well as on the origin of the 
Znr symmetry, which we do not discuss in this paper. On the other hand, the operator 
relevant to the /x-term is also written in terms of {X): 

2a— h—h+rN 
( (X) \ 2a + nN 

W ^ Y: < y MpHH, (93) 

n>-2a/N J 

where r is the minimum integer which gives 2a — h — h + rN > 0. Then, we naturally 
expect that the above operator is also dominated by the term with n = Uq: 

2a~h~h+rN 

W ^ ^n,y-^j MpHH. (94) 

Though it is possible that the fractional power which gives rise to the //-term is different 
from the one responsible for the (W), we argue that it is unnatural and assume that both 
of them are dominated by the term with same n — Uq. Then, from Eqs.(92) and (94), we 
obtain the following expression of the //-term, 

2a-h-h+rN 

W (^] " MpHH, (95) 
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which leads to Eq.(l) with the y given in (2). The same argument is also applied to the 
cases of baryon- and lepton- number violating operators discussed in the text. 
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